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Abstract 

We present an electronic model with long range interactions. Through 
the quantum inverse scattering method, integrability of the model is estab- 
lished using a one-parameter family of typical irreducible representations of 
gl{2\l). The eigenvalues of the conserved operators are derived in terms of 
the Bethe ansatz, from which the energy eigenvalues of the Hamiltonian are 
obtained. 



1 Introduction 

The Quantum Inverse Scattering Method (QISM) [|l]] is one of the most powerful 
tools in the exact study of quantum systems. It can be applied in a number of 
contexts, including both one-dimensional systems with nearest neighbour interac- 
tions such as the Heisenberg [|l]] and Hubbard [||] models, and also for the analysis 
of models with long range interactions such as the Gaudin Hamiltonians Q and 
extensions [0] . These latter constructions in particular have received renewed at- 
tention as it has been realised that the reduced BCS model, which was recently 
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proposed to describe superconducting correlations in metallic grains of nanoscale 
dimensions can be shown to be integrable through the use of Gaudin Hamil- 
tonians in non-uniform external fields Formulating the reduced BCS model 
in the framework of the QISM reproduces the exact solution originally obtained 
by Richardson and Sherman [0], and opens the way for the calculation of form 
factors and correlation functions . 

Motivated by this result, one can investigate to what extent the construction 
can be generalised to yield new classes of models, with some examples already 
given in Here, we will consider a case where an underlying superalgebraic 
structure (i.e., one with both bosonic and fermionic degrees of freedom) is em- 
ployed to yield an electronic model. The supersymmetric formulation of inte- 
grable systems can be traced back to the work of Kulish [|T0|], and recently super- 
symmetric Gaudin Hamiltonians have been analysed in detail in [[TT|]. 

In this article we present a Hamiltonian derived through the QISM from a 
solution of the Yang-Baxter equation (YBE) associated with a typical irreducible 
representation of the Lie superalgebra gl{2\l), which has the explicit form 

j j,k o-=± 

Above, the energy levels ej are two-fold degenerate, g is an arbitrary coupling 
parameter and D is the total number of distinct energy levels. Also, rij is the 
fermion number operator for energy level ej and for the parameters aj, j = 1...D 
we define 

with Xj = \Jaj/ [aj + 1). The c^, cJ., a = ±, are two-fold degenerate Fermi 

annihilation and creation operators. 

The Hamiltonian has a similar form to the reduced BCS model [Q]. There, 
Cooper pairs are scattered into vacant energy levels while the one particle states 
are blocked from scattering. In the Hamiltonian above, there is correlated scatter- 
ing depending on the occupation numbers. One of the features of this model is that 
the scattering couplings can be varied through the choice of the parameters a^. Via 
the algebraic Bethe ansatz method and using the minimal typical representation of 
gl{2\l), from which these free parameters arise, we establish the exact solvability 
of the model. Here we outline the necessary definitions and constructions, while 
full details will be presented elsewhere. 

The Lie superalgebra gl{2\l) has generators E^, i,j = 1, 2, 3 with supercom- 
mutator relations 

Above, the BBF grading [1] = [2] = 0, [3] = 1 is chosen and the elements are 
realised in terms of the Fermi operators through (cf. [[T^]) E^ = = c^c_, 
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Ef = S = cLc+, El = —a — n+, E^ = —a — n_, E^ = 2a + n, 
El = Ql, El = Ql, El = Q+ and E| = Q_, and we set = (n+ - n_)/2. 
The Casimir invariant of the algebra, C = Y.f,j=iEj Ej(—iy^\ which com- 
mutes with all the elements of gl(2\l), will also be needed, and has the eigen- 
value = —2a{a + 1) in the above representation. Below we let V{a) de- 
note the four-dimensional model on which the representation acts, with the basis 

|o). 



2 The Yang-Baxter equation and integrability 

To construct the model, we use the supersymmetric formulation of the QISM [[T0|]. 
We take the following solution of the YBE which acts onW ®W ® V{a), where 
W denotes the three-dimensional vector module of gl{2\l), 

Ruiu - f )Li3(M)L23(f ) = L23{v)Li3{u)Ri2{u - v) (2) 

with 

R{u)=I^I + ^ Y: {-lpe'::®el, (3) 

m,n=l 



and the L-operator is given by 



L(u) = /®/ + ^ ^ (-1)^6-®^;;,, (4) 



^ m,n=l 



The representations taken for the operators E^^ are as stated above, the variable u 
represents the rapidity i] is arbitrary and / is the identity operator. 

By the usual procedure of the QISM, we define a transfer matrix acting on the 
D-fold tensor product space (for distinct a^) V{ai) V{a2) ■■■ V{aD) via 

t(u) = stro {GqLqd{u - eD)--Loi{u - ei)) , 

which gives a mutually commuting family satisfying [t{u), t{v)] = 0. Above, 
strQ denotes the supertrace taken over the auxiliary space labelled by and G can 
be any matrix which satisfies [R{u), G G] =0. 

For the BBF grading we choose G = diag(exp(/5?7), exp(/5?7), 1) and by em- 
ploying the algebraic Bethe ansatz method the eigenvalues of^ the transfer matrix 
are found to be (cf. [[TOp) 

A{u) = eMf3v)I[(l--r^)l[a{v,-u) 
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^ ( mi \ ^ ^ 

+ exp(/5r/) n u - / X n ~ n '^(^-^ ~ 
-n(i-^)n<.(7.-). (5) 

where a{u) = l + r]/u. The parameters vi, Wj satisfy the Bethe ansatz equations 

TT f ^ /r«(^i-f^i) - ei - 2r]ai ( o \^ t n 

n "(T* -".)=- n ^^^;;;^y n ^, = "MPvmad. -v,). 

We now introduce the operators 

Tj = Jim ^^-^^fi^t{u), which satisfy [T,, T^] = 0. (6) 
By taking the quasi-classical expansion Tj — tj + o{rf), this leads to 



where 61 = Em,„ ® and = Ef. It is easily deduced that these 

operators satisfy [r-, , r^] = 0. 

Writing K = E«^'j('5+5'7 + S^S^ + 25/51), which satisfies \K, r^] = 0, V j, 
we define the Hamiltonian as follows; 

^ D D 1^ 1^ 1^ 

= ^ E E ^i*: - ^ E(l + 2/?e.)V',- + ^ E V^^V',- - ^ E ^ji^j + 1) 

-^ + 2Ee,K + i). 

The term involving djk may be simplified using the Casimir invariant and the 
commutation relations of the algebra gl{2\l) 

D D D D 

3 ki^3 k,j j 

= ^ - ^E^.^fc + E^. - 2E E Q%Qk. + 2j:^3i^3 + !)• 

3,k 3 3,k o-=± 3 



For g = 1/(3 we obtain the Hamiltonian (IT]), which establishes integrability since 
[H,Tj] = 0, Vj. 

From (^J^) we obtain the eigenvalues of Tj for the BBF grading, 



^ 1 ^ aa- 

A, = -2/?a, + a,5:^— -2^:^, (7) 

as the quasi-classical limit of the eigenvalues of the transfer matrix. The corre- 
sponding Bethe ansatz equations are 

Pi Ml P ^ 

/3 + E^ = E^^' E^ = 2E^- (8) 

i 7«-^j Te*-7z t^jV,-Vj 

For a given solution of the Bethe ansatz equations we find that the number of 
electrons, N = 2D — M, — n_ = M — 2P and the eigenvalue of K reads 

^K = l{M-2P){M-2P + 2). (9) 
The energy eigenvalues can be computed using (0j8|,^) and are given by 

D M D 

E = 2Y.e,-Y.^i-2gY.a,-gM. 



Similar results have been obtained for the FBB and BFB gradings, which will 
appear elsewhere. 
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